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We  present  results  of  an  experimental  study  of  the  timing  and  phase  dynamics  in  a  mode-locked  Ti:sapphire 
laser.  By  measuring  the  response  of  two  widely  spaced  comb  lines  to  a  sinusoidal  modulation  of  the  pump 
power,  we  determine  quantitatively  the  response  of  both  the  central  pulse  time  and  the  phase.  Because  of 
the  distinct  response  of  the  pulse  energy,  central  frequency,  and  gain  to  the  modulation,  we  are  able  to  dis¬ 
tinguish  their  contributions  to  the  timing  and  phase  dynamics.  ©  2007  Optical  Society  of  America 
OCIS  codes:  140.4050,  190.5530. 


The  ability  to  control  the  phase  of  femtosecond  pulses 
from  mode-locked  lasers  [1]  has  advanced  optical  me¬ 
trology  [2]  and  enabled  optical  clocks  [3].  These  ap¬ 
plications  rely  on  the  laser’s  narrow  comb  linewidth, 
whose  magnitude  depends  mainly  on  the  timing  and 
phase  jitter  of  the  pulse  train.  Recently,  there  has 
been  interest  in  how  the  dynamics  of  mode-locked  la¬ 
sers  affect  the  properties  of  the  comb  [4-6].  Cur¬ 
rently,  technical  noise  is  the  main  source  of  linewidth 
broadening,  but  the  fundamental  limit  is  set  by  am¬ 
plified  spontaneous  emission  [7].  The  quantum-noise- 
limited  comb  line  shape  has  not  been  determined. 
While  the  statistical  properties  of  the  quantum-noise 
sources  are  well  understood,  the  linear  response  of 
the  mode-locked  pulse  parameters  to  perturbations 
must  also  be  known  to  calculate  the  line  shape.  This 
linear  response  has  not  been  known,  except  qualita¬ 
tively,  and  in  this  paper  and  an  earlier  paper  [8]  we 
directly  measure  this  response.  In  our  earlier  work, 
we  showed  that  the  gain  must  be  included  among  the 
parameters  whose  response  is  considered,  and  we 
quantitatively  measured  the  response  of  the  pulse 
energy,  the  central  frequency,  and  the  round-trip  gain 
to  a  perturbation.  Here,  we  measure  the  coefficients 
that  couple  those  three  parameters  to  the  central 
pulse  time  and  the  phase.  With  these  results,  we 
have  enough  information  to  fully  characterize  the  re¬ 
sponse  of  the  laser  pulse  parameters  to  perturbations 
and  hence  to  predict  the  comb  line  shape  of  the  laser 
studied.  In  addition,  this  linear  response  plays  a  key 
role  in  determining  the  pulse  stability.  Thus,  this 
study  also  has  important  implications  for  laser  de¬ 
sign. 

Previously,  we  used  a  sudden  modulation  of  the 
pump  power  and  studied  the  response  of  the  pulse 
energy  and  central  frequency  in  the  time  domain  [8]. 
Here,  we  use  a  sinusoidal  modulation  and  examine 
the  response  of  the  laser  parameters  in  the  frequency 
domain.  We  have  verified  that  such  a  swept  sine  tech¬ 
nique  gives  the  same  results  for  the  energy  and  cen¬ 
tral  frequency  as  in  [8] .  Because  the  sinusoidal  wave¬ 
form  has  no  low-frequency  components,  it  does  not 


affect  the  feedback  loop  central  to  the  technique  de¬ 
scribed  here  for  measurement  of  timing  and  phase. 
When  the  pump  power  of  the  laser  is  modulated  at 
frequencies  near  the  relaxation  oscillation  frequency 
(approximately  400  kHz  in  the  laser  studied  here), 
the  gain,  intensity,  and  central  frequency  all  have  dif¬ 
ferent  responses.  Thus,  we  may  distinguish  their  con¬ 
tributions  to  changes  in  the  central  pulse  time  and 
phase. 

The  pulses  are  characterized  by  four  parameters: 
the  pulse  energy  w  =  weq  +  Aw,  the  central  frequency 
in  =  TTTeq+ Atjt,  the  central  pulse  time  r=Teq+Ar,  and 
the  phase  0=0eq+A0,  where  xeq  denotes  the  equilib¬ 
rium  value  of  each  quantity  and  Ax  its  change.  Be¬ 
cause  of  the  system’s  invariance  under  time  and 
phase  translations,  we  may  choose  t6„  =  0.  We  choose 
0eq=2TrfoT,  where  f0  is  the  equilibrium  carrier- 
envelope  offset  frequency.  Following  common  practice 
in  soliton  theory,  the  phase  is  defined  with  respect  to 
a  common  reference  rather  than  the  pulse  envelope. 
With  this  definition,  movement  of  a  mirror  in  the  la¬ 
ser  cavity  results  in  changes  in  both  r  and  9,  and  the 
carrier-envelope  phase  slip  A0ce  is  related  to  Ad  and 
At  by  Ad=-Adce-MnAr.  To  these  four  pulse  param¬ 
eters  we  add  the  round-trip  gain  g=g(,q  +  Ag  as  a  dy¬ 
namic  variable  [8] .  The  linear  response  of  the  laser  to 
perturbations  is  then  [9] 

dv 

—  =  -A-v  +  S,  (1) 


where  v  =  (Ag ,  Aw ,  Am ,  At,  A0)( ,  S  is  a  vector  of  noise 
sources  for  each  parameter,  and  A  is  a  matrix  of  co¬ 
efficients  that  describe  the  linear  response  of  each 
parameter  to  changes  in  itself  or  the  others.  Once  A 
and  S  are  known,  one  can  calculate  the  timing  and 
phase  jitter  [9] ,  which  leads  to  the  optical  frequency 
comb  linewidth  [10,11].  Here,  we  measure  the  ele¬ 
ments  of  A  that  couple  w,  g,  and  m  to  t  and  0.  As  t 
and  9  cannot  drive  w,g,  and  tct  due  to  time  and  phase 
invariance  of  the  system,  AXT  and  Axe  are  zero. 
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The  nonthermal  change  in  the  repetition  rate  /'rep 
caused  by  a  few  percent  modulation  of  the  pump 
power  is  less  than  1  Hz,  making  an  electronic  mea¬ 
surement  of  the  pulse  timing  dynamics  over  a  few  {is 
extremely  difficult.  Microstructured  fiber  used  for 
self-referenced  detection  [1]  of  the  carrier-envelope 
offset  frequency  f0  converts  intensity  fluctuations 
into  phase  (and  thus  frequency)  fluctuations  [12], 
complicating  the  use  of  that  technique  for  measure¬ 
ment  of  the  phase  dynamics.  To  circumvent  both  of 
these  problems,  we  instead  measured  the  response  of 
comb  lines  by  monitoring  the  heterodyne  beat  be¬ 
tween  the  laser  being  perturbed  and  a  reference  la¬ 
ser.  The  multiplication  up  to  optical  frequencies  mag¬ 
nifies  small  changes  in  t.  By  measuring  the  response 
of  two  widely  separated  comb  lines,  we  can  separate 
the  contributions  from  changes  in  t  and  6. 

The  linear  cavity,  Kerr  lens  model-locked  Ti:  sap¬ 
phire  laser  under  study  had  a  repetition  rate  of 

93.5  MHz,  with  dispersion  compensation  provided  by 
double  chirped  mirrors  and  CaF2  prisms.  The  spec¬ 
trum  was  centered  at  826  nm  with  a  width  of  80  nm. 
It  was  pumped  with  a  532  nm  diode-pumped  solid 
state  laser  that  passed  through  an  acousto-optic 
modulator  that  could  deflect  up  to  5%  of  the  power 
with  a  bandwidth  of  approximately  5  MHz.  The 
pump  power  before  the  acousto-optic  modulator  was 

5.5  W.  The  modulation  frequency  was  swept  from 
50  to  800  kHz  and  the  modulation  depth  was  ap¬ 
proximately  0.3%.  We  measured  the  pump  power  as  a 
function  of  time  to  normalize  by  the  modulation  am¬ 
plitude  and  phase.  We  also  captured  the  laser  power, 
spectrally  resolved  with  a  monochromator,  as  a  func¬ 
tion  of  time  using  a  photodiode.  The  integrated  spec¬ 
trum  and  centroid  give  the  pulse  energy  and  central 
frequency,  respectively. 

We  locked  the  tenth  harmonic  of  the  repetition  rate 
of  the  laser  under  study  to  the  935  MHz  repetition 
rate  of  a  Ti:  sapphire  ring  laser.  The  bandwidth  of  the 
feedback  loop  was  restricted  to  10  kHz  to  ensure  that 
the  loop  could  not  compensate  for  changes  in  fre p  from 
the  modulation  in  the  pump  power.  The  actuator  for 
the  cavity  length  was  a  piezo  mounted  on  the  high  re¬ 
flector  of  the  laser  under  study.  Light  from  both  la¬ 
sers  was  coupled  into  single  mode  fibers  and  com¬ 
bined  to  generate  an  optical  heterodyne  beat.  For  two 
unperturbed  lasers  with  harmonically  related  repeti¬ 
tion  rates,  the  beat  frequency  is  the  difference  in  f0 
(modulo  the  lowest  repetition  rate)  between  the  two 
lasers.  The  carrier-envelope  offset  frequencies  of  both 
lasers  were  free  to  drift  but  did  not  change  signifi¬ 
cantly  over  a  few  {is,  the  time  scale  of  the  experi¬ 
ment.  We  captured  the  heterodyne  beat  with  an  os¬ 
cilloscope  and  used  a  nonlinear  fitting  routine  to  find 
the  beat  frequency  as  a  function  of  time. 

The  frequency  of  the  nth  comb  line  is  /„  =/i/n!p+/0. 
The  change  in  the  frequency  of  the  heterodyne  beat  is 
equal  to  the  absolute  value  of  the  change  in  the  comb 
line  frequency  |A fn\.  A  change  in  r  leaves  the  central 
comb  line  near  i3Teq  fixed  and  moves  the  nth  line  by  a 
value  proportional  to  -/0,  where  fc  is  the  central 
frequency,  assumed  here  to  be  the  centroid  of  the 
spectrum.  A  change  in  6  affects  all  comb  lines  equally. 


Given  changes  in  two  comb  lines  of  frequency  fn  and 
fm,  the  rate  of  change  of  t  is 


Af= 


A fn  -  kfrn 


(2) 


The  further  apart  the  comb  lines,  the  more  sensitive 
the  measurement  of  r.  The  rate  of  change  of  6  is 

bfnifm  -fc-  fo)  -  bfjfn  ~fe~  fo) 

A0=2tt - .  (3) 


We  used  a  monochromator  with  a  bandwidth  of  1  nm 
to  select  two  groups  of  comb  lines  at  350  and 
390  THz.  There  is  an  ambiguity  in  the  sign  of  the 
change  in  fn  and  fm  because  of  the  beating  down  to 
radio  frequencies.  Results  are  shown  in  Fig.  1. 


Fig.  1.  Experimental  results  for  (a)  pulse  energy  w,  (b) 
central  frequency  tu,  (c)  central  pulse  time  r,  and  (d)  phase 
6.  The  amplitude  (triangles)  and  phase  (squares)  of  the  re¬ 
sponse  of  each  variable  are  plotted,  as  well  as  fits  (solid 
curves)  in  (c)  and  (d)  to  Eq.  (4)  using  the  measured  re¬ 
sponse  of  w  and  tu. 
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Under  the  assumption  that  changes  in  0  and  t  do 
not  affect  the  other  variables,  their  response  is 

A  HU) 

— —  =  -Anuw o(U)  -A^g0(n) -Almoin), 

J-R 

A  6»(H) 

— —  =  -A,)wW0([l)  -A^0(n) -Atem0(n),  (4) 

J-R 

where  w0,  g0,  and  ur0  are  the  complex  amplitudes  of 
those  variables  as  functions  of  the  angular  modula¬ 
tion  frequency  11.  Therefore,  one  can  extract  the  Axy 
coefficients  by  fitting  the  measured  r0  and  0{)  to  w0, 
g0,  and  thq.  We  used  the  measured  w0  and  m0,  shown 
in  Figs.  1(a)  and  1(b),  and  used  a  theoretical  model 
for  g0  (the  frequency  domain  version  of  what  was  de¬ 
scribed  in  [8]). 

We  assigned  the  sign  of  the  extracted  Axy  coeffi¬ 
cients  using  the  known  sign  of  the  dispersion  using 
ATm=-D/TR,  where  D  =  -(3"  is  the  average  dispersion 
in  the  cavity  [9],  We  find  ATCT  =  -4.0X  1010  as/(s  THz), 
Ara)  =  2.0x  109  as/(s  nJ),  and  Ato=9.4X  105  rad/ 
(s  nJ).  The  data  are  consistent  with  the  gain’s  having 
no  direct  effect  on  either  r  or  9  and  with  the  central 
frequency  having  no  effect  on  9,  i.e.,  Affm=Agg=A7g 
=  0.  The  coefficients  coupling  pulse  energy  to  r  and  9 
are  related  to  the  Kerr  nonlinearity  in  the  Ti:  sap¬ 
phire  crystal.  The  Kerr  phase  shift  leads  to  a  nonzero 
Ato,  while  Kerr  shock  [13,14]  contributes  to  a  non¬ 
zero  Am.  Other  processes  may  also  contribute. 

The  variation  of  f0  with  pump  power  is  often  used 
to  lock  it  [15,16].  However,  many  competing  pro¬ 
cesses  are  involved,  some  of  which  affect  t  and  others 
6  (both  play  a  role  in  the  phase  evolution  because 
A0ce=-A6+TzAT).  In  addition  to  the  direct  changes  in 
phase  and  timing  caused  by  the  Kerr  effect  (the  Kerr 
phase  shift  and  shock,  respectively),  gain-  and 
energy-driven  frequency  pulling  changes  m  [8] ,  which 
changes  r  through  dispersion.  The  technique  de¬ 
scribed  here  allows  quantitative  measurement  of 
these  different  processes.  Using  d9ce/dw=ASw-mAnv, 
we  note  that  we  observe  a  partial  cancellation  of  the 
intensity  related  effects  to  A fn,  as  predicted  previ¬ 
ously  [13,14]. 

In  summary,  we  have  developed  a  comb-based  tech¬ 
nique  for  sensitive  measurement  of  the  linear  re¬ 
sponse  of  the  timing  and  phase.  By  taking  advantage 
of  the  differing  response  functions  of  the  pulse  energy 
and  central  frequency,  we  were  able  to  fully  charac¬ 
terize  the  timing  and  phase  response  of  the  laser. 
These  results,  along  with  those  presented  in  [8],  pro¬ 


vide  a  quantitative  measurement  of  the  linear  re¬ 
sponse  matrix  A.  An  accurate  calculation  of  the 
quantum-noise-limited  comb  line  shape  for  this  laser 
is  now  possible.  Additionally,  we  expect  that  the 
unique  ability  of  the  technique  described  here  to  dis¬ 
tinguish  between  changes  in  the  timing  and  phase 
and  identify  the  sources  of  these  changes  will  prove 
useful  in  the  study  of  other  types  of  mode-locked  la¬ 
sers. 
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